Diffusive transport in graphene: the role of interband correlation 
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We present a kinetic equation approach to investigate dc transport properties of graphene in the 
diffusive regime considering long-range electron-impurity scattering. In our study, the effects of 
interband correlation (or polarization) on conductivity are taken into account. We find that the 
conductivity contains not only the usual term inversely proportional to impurity density Ni, but 
also an anomalous term that is linear in Ni. This leads to a minimum in the density dependence 
of conductivity when the electron density N c is equal to a critical value, N c . For iV c > N c the 
conductivity varies almost linearly with the electron density, while it is approximately inversely 
proportional to N c when N c < N c in the diffusive regime. The effects of various scattering potentials 
on the conductivity minimum are also analyzed. Using typical experimental parameters, we find 
that for RPA screened electron- impurity scattering the minimum conductivity is about 5.1e 2 /h 
when 7V C « 0.32^. 

PACS numbers: 81.05.Uw, 72.10.Bg, 73.40.-c 



I. INTRODUCTION 

Recently, graphene has attracted a great deal of ex- 
perimental and theoretical interest ±&2d£JL In this two 
dimensional system, low energy electrons behave as mass- 
less relativistic fermions due to their linear energy spec- 
trum around two nodal points in the Brillouin zoneJ 
Such unusual electronic properties lead to high mobility 
as well as a long mean free path at room temperature, 
making graphene a promising candidate for future elec- 
tronic applications. 8 

In the first experiments involving graphene, it was 
found that there exists a finite "residual" conductivity, 
er res , in the carrier-density dependence of conductivity 
at zero gate voltage, and its value is about 4e 2 /h. Fur- 
thermore, it was found that the conductivity varies lin- 
early with carrier density when it is large. Much theoret- 
ical effort has been devoted to quantitatively explain the 
observed "residual" longitudinal conductivity. Actually, 
the existence of such a " residual" conductivity in perfect 
(scattering-free) single-layer graphene was predicted long 
before its experimental confirmation ! 9 ' 10 ! 11 ! 12 Ludwig et 
al. obtained different minimum values, cr rcs , using two 
different approaches: a xes — (ir/2)e 2 /h using the Kubo 
formula and cr rcs = (4/7r)e 2 //i using a definition of con- 
ductivity in the sigma models In many recent works, 
the issue has been further confused with findings of even 
more values of residual conductivity. In Ref s . IT3lfl4lfl~5lfl6l 
er res = (4/7r)e 2 //i was obtained, while Ziegler predicted 



ne 2 /h 
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(n/2)e /h was also obtained in 



Ref.[18|. More recently, Ziegler demonstrated that all 
these values of a Tes can be obtained within the Kubo for- 
mulism by taking different orders of the zero- frequency dc 
and zero dissipation limits. 19 Performing numerical cal- 
culations with the Kubo formula, Nomura and MacDon- 
ald obtained cr rcs = (4/7r)e 2 //i in the case of short-range 
scattering and cr xes = 4e 2 /h for the Coulomb scattering 



case 
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Employing the Boltzmann equation, Adam, et al. 
demonstrated that the minimum "residual" conductiv- 
ity arises from nonvanishing electron density at zero gate 
voltage which may be induced by impurity potentials. 21 
Analyzing random fluctuations of gate voltage, they 
found that the value of cr rcs is not universal but depends 
on the impurity concentration. Apart from this partic- 
ular issue, the observed almost-linear variation of con- 
ductivity with electron density can be easily understood 
within the Kubo formula framework 2 ^ as well as with the 
Boltzmann equation. 21 

In this paper, we present a kinetic equation approach 
to investigate transport in graphene considering long- 
range electron(or hole)-impurity scatterings. Here, in- 
terband correlations (polarization effects) are taken into 
account, whereas in all previous studies, such interband 
correlations associated with electron/hole-impurity scat- 
terings were ignored. We find that the conductivity in 
graphene contains two terms: one of which is inversely 
proportional to impurity density, while the other one 
varies linearly with the impurity density. This results 
in a minimum (rather than " residual" ) conductivity at a 
nonvanishing critical electron density, N c , in the electron- 
density dependence of conductivity. For electron density 
N e larger than N c , the conductivity increases almost lin- 
early with increasing N c , while, the conductivity is ap- 
proximately inversely proportional to N e for N < N c . 
We also demonstrate the effects of various scattering 
potentials on the conductivity minimum. Considering 
RPA screened electron- impurity scattering, we find that, 
for typical experimental parameters, the critical electron 
density is about 0.32A^ (Ni is the impurity density) and 
the value of minimum conductivity is equal to 5.1 e 2 /h. 

The paper is organized as follows. In Sec. II the ki- 
netic equation for noncquilibrium distribution functions 
as well as its solution are presented. Also, the conductiv- 
ity is exhibited in terms of microscopically derived relax- 
ation times. In Sec. Ill we present our analytical results 
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for the conductivity for several different scattering po- 
tentials. Finally, the conclusions are summarized in Sec. 
IV. 



II. KINETIC EQUATION AND SOLUTION 

A. Kinetic equation 

In the Brillouin zone of graphene, there are six 
points at which the energy of carriers vanishes and 
the conductance band touches the valence band: p = 
(±4 7 r/(3v / 3a),0), p = (27r/(3V3a), ±27r/(3a)), and p = 
(— 27r/(3v3a),±27r/(3a)) with a as the lattice spacing. 
These points correspond to two inequivalent Dirac nodes, 
K and K' . In present paper, we are interested in 
the transport of carriers in graphene with momenta 
near these Dirac points. The Hamiltonian of an elec- 
tron with two-dimensional momentum, p = (p x ,Py) = 
(pcos0p,psin0 p ), near the v = K or K' Dirac nodes 
can be written as 



(p) = ll&xPx + sgn(v)a y p v ] 



(1) 



with sgn(u) = 1 or —1 for v = K or K' , and 7 = 
V3aa/2 is a material constant (a is the hopping param- 
eter in tight-binding approximation). 7 is equal to the 
Fermi velocity, which is independent of carrier density in 
graphene. 

The Hamiltonian ([1]) can be diagonalized, resulting in 
two eigen wavef unctions, </?^p(r) = it M (p)e lp ' r , and two 
eigenvalues, e M = (— l) M+1 e p , with /i = 1, 2 as the helicity 
index, e p = jp, and 



>) 



(P) 



1 

V2 



,(-1) 



l) M isgn(u)0 p 
1 
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£i(p) and £2(p) are just the dispersion relations of the 
conduction and valence bands, respectively. 

It is useful to introduce a unitary transformation, 
Up = (ui(p), 1*2 (p)), which corresponds to a change from 
a pseudospin basis to a pseudo-helicity basis. Applying 
this transformation, the Hamiltonian (TTJ) is diagonalized 



as h { v) (p) = [U p v ']+h^(p)U p v > = diag[ £l (p),e 2 (p)]. 
Note that Kq (p) is independent of node index v, in- 
dicating the existence of a valley degeneracy. 

In a realistic graphene system, the carriers experi- 
ence scattering by impurities. We assume that, in the 
pseudospin basis, the interaction between carriers and 
impurities can be characterized by an isotropic poten- 
tial, V(|p — k|), which corresponds to scattering a car- 
rier from state p to state k. In the pseudo-helicity ba- 
sis, the scattering potential takes the transformed form, 

We are interested in the current in a graphene system 
driven by a dc electric field, E. In the pseudo-spin basis, 
this electric field can be described by a scalar potential, 
V = eE ■ r, with r as the carrier coordinate. From Eq. ((T|) 



it follows that the pseudo-spin single-particle current op- 
erator, j^ v \p), has vanishing diagonal elements: 



jW( P ) = v P C(p) 



(v), 



(3) 



The observed net current, given by J = 
e^E w ,p^Ej W (p)P M (p)] _ with pM(p) as the 
pseudospin-basis distribution function and g s as 
the spin degeneracy of graphene, can be also determined 
in pseudo-helicity basis via 



J = 9s e^Tr[jW( P )^)(p)], 



(4) 



with j W (P) 



[U p v '] + p^ (p)Up^ being the pseudo-helicity-basis single- 
particle current operator and distribution function, re- 
spectively. Explicitly, Eq. (0| can be rewritten as 



and p^Xp) 
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+2sgn(v)[pxn]Im[[^>]i 2 (p)]} 



(5) 



with [p^]ij,u(p) as the elements of the distribution func- 
tion. In the derivation of Eq. ([5|) , the Hermitian feature 
of the distribution function, i.e. p( v \p) = [p^] + (p), has 
been used. From Eq. it is evident that contributions 
to current arise not only from the diagonal elements of 
the distribution function, but also involve its off-diagonal 
elements. 

To carry out the calculation of current in graphene, it 
is necessary to determine the carrier distribution func- 
tion. In the pseudo-spin basis, the kinetic equation for 
the distribution, p^(p), can be written as 

eE-V p p^(p)+i[h^(p) ) p^(p)} = -/», (6) 
with /w as the collision term. Applying the unitary 

(v) 

transformation U p , the kinetic equation for the dis- 
tribution in the pseudo-helicity basis, p^(p), takes the 
form 



eE-<{VppM(p) 



p(p),v p (c/w) u& 
i[h£ ) ( P ),pW{p)] = -iw, (7) 



with /(») = [U p v) }+I^[U p v) ]. 

It should be noted that the second term on the left- 
hand side of Eq. ([7]) is associated with interband tun- 
neling since the off-diagonal elements of matrix eE • 



(v) 



p are just the interband-tunneling matrix 

elements < (p)|eE • r\u^ (p) > (p = 3 — p) with r as 
carrier coordinate i 23 ' This term of Eq. ([7]) results in a 
component of [/3^]i2(p) which depends on the strength 
of electric field via e~ b / E (b is a real, positive parame- 
ter), and is much less than unity in the linear response 
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regime. Thus, the interband tunneling term of Eq. |7J) is 
not actually linear (despite its formal appearance), and, 
in fact, it can be ignored at low fields. Hence, Eq. |0 can 
be rewritten as 

eE- V p/ 5W(p) + z[^' ) (p),^)(p)] = (8) 
with l( v > as the collision term given by 

f(v) _ yr p< , y< pa _ pr y< _ p< ya /q\ 
vp vp ' vp up vp vp vp up' \ J 

Here, G%£' < and respectively, are the nonequilib- 

rium Green's functions and self-energies for carriers near 
node v = K or K' . 

In the kinetic equation above, electron-impurity scat- 
tering is embedded in the self-energies, £™ ,<: - In the 
present paper, we only consider electron-impurity col- 
lisions in the self-consistent Born approximation. It is 
widely accepted that this is sufficiently accurate to an- 
alyze transport properties in the diffusive regime.— Ac- 
cordingly, the self-energies take the forms: 



^vp 



k 



T„( P ,k)GX^T+( P ,k) 



(10) 



In the present paper, we restrict our considerations to 
the linear response regime. In connection with this, all 
the functions, such as the nonequilibrium Green's func- 
tions, self-energies, and distribution function, can be ex- 
pressed as sums of two terms: A = Aq + Ax, with A 
representing the Green's functions, self-energies or dis- 
tribution function. Aq and Ax, respectively, are the un- 
perturbed part and the linear electric field part of A. In 
these terms, the kinetic equation for the linear electric 
field part of the distribution, (p) , takes the form 

e E.V p ^(p) + «[^,pW(p)]^/W, (11) 

with Iv as the linear electric field part of the collision 
term I v . This equation can be further rewritten explicitly 

as 



eE.V p ^U(p) = -[/«] 



and 



(«)- 



12 (,P) 



■[I 



(1) 



(12) 



(13) 



with (px)nu(p) and (po)nn(p), respectively, as the ele- 
ments of px(p) and Pq (p) = diag[n F (ei(p)), np(ei(p))], 
and n is a unit vector perpendicular to the graphene 
plane. 

To further simplify Eq. (jlip . we employ a two-band 
generalized Kadanoff-Baym ansatz (GKBA) * 25 ' 26 This 
ansatz, which expresses the lesser Green's function in 
terms of the Wigner distribution function, has been 
proven sufficiently accurate to analyze transport and op- 
tical properties in semiconductors. 27 To first order in the 
dc field strength, the GKBA reads, 

G<(p,c) = -G5(p,w)ft(p) +Pi(p)G§(p,w), (14) 



where the unperturbed retarded and advanced Green's 
functions are diagonal: GQ ,a (p,w) = diag[(w — £i(p) ± 
i<5) _1 ,(a> — £2(p) ± iS)^ 1 ]. In our treatment, the ef- 
fect of G[' a (p,uj) on the distribution function has been 
ignored because these linear electric field parts of the 
retarded and advanced Green's functions lead to a col- 
lisional broadening effect on px{p), which plays a sec- 
ondary role in transport studies. Further, in our treat- 
ment, we ignore intervalley transitions of carriers between 
different Dirac points since the corresponding rates are 
very small in the presence of carrier-impurity scattering. 

Within the framework of these considerations, the scat- 
tering term Iv in Eq. (fTTj) can be expressed in terms of 

the distribution function: its diagonal elements, [J^]^, 
can be written as 



f/ (1) l - ttN- 



k 



|y( P -k)| 2 %(")(p)- e W(k)] 



x {[1 + cos^ p - M]{[p[ v) U(p) ~ [p^Wtk)} 

-(-l)"sgQ(«) sin(0 p - k ) lm{[p{ v) } 12 (p) + [p^] 12 (k)}}l5) 
while the off-diagonal element, [!^]i2, takes the form 

= ^Y. I"<P - k)l J «[4"»(P) - 4"'(k)] 

x{( 



-l)" +1 iBgn(w)siii(^p-^ k ){^ 

[1 - cos(0 p - M] {[p^Mp) ~ ^] 2 i(k)}}(16) 



In the derivation of these equations, the effects of the real 
parts of the retarded and advanced Green's functions on 
Iv 1 ' 1 have been ignored. 

B. Conductivity and solution of the kinetic 
equation 

Since e Al (p), as well as the equilibrium distribution 
Po (p) > depend only on the magnitude of momentum, 
the dependence of p^ (p) on momentum angle can be 
evaluated explicitly. From Eq. (fTTj) we see that the di- 
agonal elements of px (p) can be written as [vj, (p) = 

Vp£< o) (p) - (-ir +1 7P/p] 



|p£ o) W(p) = eE.vW(p)A£0(p), 



(17) 



while the off-diagonal element, [pi ] n (p) , takes the form 



[/5 i ' ;) ]i 2 (p) = ^[Exp.n]$W(p). 



V 



(18) 



The functions A^(p) and $^(p) depend only on the 
magnitude of momentum and are determined by coupled 
equations: 

de), (p) >(p) L > 



4 



{sgn(v)A^(p)-lm[^ v \p)]}, (20) 



and 



2 7P Im[$(' u )(p)] = J2 -i W W (p)]| , (21) 

with t^ (p) as microscopically determined relaxation 

times independent of node index v: b"/^ '^(p)] 1 = 
^E k l^(p-k)| 2 5[ £A1 (p)- £Al (k)]A( Q . b )(0 k );AW(0) = 
sin 2 <f>, and A^{<j>) = (1 — cos0) 2 . Solving the algebraic 
equations (fl"9)) - ([2"Tj) yields explicit analytic expressions for 

the functions and & v \p) in terms of the relaxation 
times as 



Re[$M(p)] 



Im[$(")(p)] 



-f^E{7 / {[4 l ] r (p)} l^) 



a4 u) (p) 



sgn(«) 

16 7 3 p 3 

*E 



rfb) r< 6 >(p) 



d {\g\M) 
del v \p) 



IP 



(23) 



and 



A< u) (p) 



W ^ } (p) 

+sgn(v)Im[$ ( " ) (p)]. 



(24) 



Substituting Eqs. (JX7J) and (JTHJ) into Eq. ©, J can be 
further rewritten as 



J = ^. 9s7 2 EE{[ A i t ' ) b)+ A 2 t ' ) b)]+2sgn(«)Im[$W(p)]} 



(25) 

Inserting the explicit forms of A^, (p) and $(^(p), we 
find that the conductivity, cr, takes the form 



E { 



4 7 2 pM 6 )(p) 
5e^( P ) 



(26) 



with [t^ (p)] 



i w (p)]- 



(P)]" 



From Eq. (j2"B"|) it is clear that, the conductivity in 
graphene involves not only a term proportional N^~ , 
but also a term linear in impurity density. How- 
ever, there is no term independent of impurity scat- 
tering. This differs significantly from previous results. 



Adams et al. and Galitski et al. found the conduc- 
tivity to always be proportional to N^ 1 f 21 i 22 while the 
other previous calculations indicated that the conduc- 
tivity always contains a term independent of impurity 
scattering^a^lS^lLia^a 

It should be noted that the anomalous conductivity 
term proportional to Ni arises from the nonvanishing of 
]i2(p)- This can be seen from Eq. ([25]) : the conduc- 
tivity depends on the imaginary part of [/5^]i2(p), which 
is proportional to [t^ (p)}^ 1 ■ Since [Pi ]ia(p) describes 
interband coherence, the anomalous conductivity term is 
a result of interband correlation. 

From Eqs. (|20p and (12"T]) it follows that the imaginary 

part of [pi ]i2(p) is proportional to Ni, while its real part 
is independent of impurity density. This can be under- 
stood from the fact that both the driving electric field and 
impurity scattering can cause transitions between two 
pseudo-hclicity states: they result in a change of carrier 
momentum but retain the pseudospin unchanged, leading 
to a change of pseudo-helicity. Obviously, the probabil- 
ity of this transition is proportional to the strength of the 
momentum change, which is determined by the electric 
field and/or by impurity scattering. Thus, impurity scat- 
tering results in a term of [pi ]i2(p) proportional to Ni, 
while the A^-independent term of [pi^i^Cp) is a result 
of the momentum change induced directly by the electric 
field. 



III. RESULTS AND DISCUSSION 

As pointed out by Adam, et at, "intrinsic graphene" 
with chemical potential (or Fermi energy) precisely at the 
Dirac point is experimentally unrealizable since charged 
impurity disorder or spatial inhomogeneity in the sys- 
tem can lead to a nonvanishing induced graphene carrier 
density^ Accordingly, all experimental graphene sam- 
ples are extrinsic since some free carriers always exist in 
the system. In connection with this, the conductivity of 
graphene arises mainly from the contribution of one type 
of carrier, electrons or holes, depending on the experi- 
mental conditions. 

In the present paper, we focus on electron transport in 
graphene at zero temperature. In this case, the chemical 
potential (or, equivalently the Fermi energy, Ep) is given 
by Ep = 7 fc.F with the 2D Fermi wavevector kp depend- 
ing on the electron density through fop = y^AirNe/ g s g v 
( we choose the spin and valley degeneracies, respec- 
tively, as g s — g v = 2 in the present paper). The zero- 
temperature conductivity is given by 



lr=o 



e —. — lk F 
4tt 



» 



(k F ) 



1 



4 7 2 fc 2 TW(£: F ) 



.(27) 



To analyze the effects of scattering potentials on cr, 
we consider short-range (SR), Thomas-Fermi (TF), as 
well as random-phase-approximation (RPA) screened 
Coulomb interactions. 
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background static dielectric constant 
3.5 3 2.5 2 1.5 1 



tential, o~\j^_ Q , reduces to 




FIG. 1: The dependencies of minimum conductivity and crit- 
ical density on r 3 for TF and RPA screened Coulomb poten- 
tials. 7 is chosen to be 1.1 x 10 J m/s. 



Short-range screened Coulomb interaction- For a short- 
range screened interaction potential, V"(p — k) be- 
comes independent of electron momentum: V(p — k) w 
7T7/ (2fcp).— In connection with this, the relaxation times 
r[ a \k F ) and r^ b \k F ), respectively, reduce to r[ a \k F ) — 

WNe/iNiTT-ykp) and T^(k F ) = 8N c /{3irNijk F ). Thus, 
the zero-temperature conductivity for a short-range po- 



iSR 
\T=0 



e 

7T 



16iV ZwNj 



irNi 32N r 



(28) 



Thomas-Fermi screened potential- We also evaluate the 
conductivity for a 2D Thomas-Fermi screened potential, 
V TF (q), given by 



2ire n{q + q s ) 



(29) 



with q = p— k, q s — 4e 2 k F / '( K 7) an d K as the static back- 
ground dielectric constant. In this case, the scattering 
times 4 a) {k F ) and T (b \k F ) can be obtained analytically 
as: 



1 



r< o) (A*)J 



TF 



Ni 
N 



■jk F G(2r s ), 



and 



W{k F ) 



TF 



N 



jk F F(2r s ), 



(30) 



(31) 



with r s = e 2 / (Atteqk,^). The functions G(x) and F{x) 
are defined as 



G(x) 



and 



X 



dO- 



(sin | + x) 

3a; — 
3a; — 



Id > 1 



^-x 2 + |x| (3x 2 - 2) arccos(l/x) [x 2 - l]" 1 / 2 

^fx 2 + x(3x 2 - 2)Re [arctanh(l/Vl - x 2 )] [1 - x 2 ]- 1 ' 2 < x < 1 ' 



-^(1-cos.) 2 



(sin ■ 



| - 4a; + 3x 2 tt - + |x 3 |(8 - 6a; 2 ) arccos(l/x)(x 2 - 1)~ 3 / 2 |x| > 1 

| - 4a; + 3x 2 tt - - x 3 (8 - 6x 2 )Re [arctanh(l/Vl - a; 2 )] (1 - a; 2 )" 3 / 2 < x < 1 



(32) 



(33) 



Thus, the conductivity for the TF potential, cr| T=0 , is 
determined by 



,TF 

It=o 



iV,G(2r s ) 



NF(2r s 
4N C 



(34) 



RPA- screened Coulomb potential- We also examine the 
conductivity of graphcne in the presence of RPA-screened 
electron-impurity scattering, with the potential taking 



the form: 



V(q) 



2-Keonqe(q) ' 



(35) 



and the dielectric function e(q) for the massless Dirac 
energy spectrum is given b y 28 ' 29 

£{q) = l+ *±\e(2k F -q) (l-$- 



+0{q - 2k F ) 



6 



4k, 



arcsin(2fci?/q r ) 



(36) 



Substituting the RPA-screened potential into the expres- 
sions for the scattering times, we obtain 



1 



r[ a \k F ) 



J RPA 



Ni 
No 



jk F G[4r s /(2-nr s )}, (37) 



and 



r( fc )(fc F ) 



RPA 



= ^ 7 k F F[4r s /(2-wr 3 )}. (38) 



Hence, from Eq. (f2"T)) . it follows that the conductivity for 
the RPA-screened Coulomb potential case is given by 



I RPA 

\t=o 



+ 



A,G[4r s /(2 - 7rr s )] 
N i F[4r s /{2--Kr a )] 



4N r 



(39) 



From Eqs.([2"8|). (f34|) . and(39|) it is evident that a 
minimum exists in the electron-density dependence of 
graphene conductivity for short-range, Thomas-Fermi as 
well as RPA-screened Coulomb potentials. The corre- 
sponding critical values of electron density, N c , and the 
minimum conductivities, er m j n , are shown in table I. Fur- 
thermore, we verify that, for N c > N c , the conductivity 
varies almost linearly with the electron density, while it 
is inversely proportional to N c when N c < N c . 

From Eq. (f28|) we see that, for SR screened scattering, 
N c /Ni and <7 m ; n are constants independent of r s . How- 
ever, for the TF and RPA potentials, they are functions 
of r s , which is, in principle, a tunable parameter through 
its dependence on the background dielectric constant k. 
In Fig. 1, we plot the dependencies of iV c and <7 m i n on r s . 
This figure enables us to determine the values of <7 m in 
and N c for various experimental samples. 

Experimentally, graphene is usually fabricated using 
a Si02 substrate and hence the background static di- 
electric constant,, k, can be estimated as k = 2. 45. 21 
The measured Fermi velocity 7 = 1.1 x 10 5 m/s leads to 
r s 0.813.— For this r s value, we obtain 



(T 

a 



TF 
min 
RPA 
min 



4A2e 2 /h, N t 



TF 



O.llJVi, 



5.lle 2 /h, iV c RPA « 0.327V;. 



(40) 
(41) 



The experimentally observed residual conductivity is 
about 4e 2 jh and our obtained minimum conductivity is 
in good quantitative agreement with it. The contention 
of Adam, et al, about finite electron density at zero gate 
voltage is plausible in describing the observed plateau 
near zero gate voltage.— However, the vagaries of ear- 
lier theories described in the Introduction, and our own 
prediction of an unlimited increase of conductivity with 
decreasing N c for N c < N c (which tends to infinity as 
N c — ► 0) calls for a re-evaluation of the assumption of 
linearity. We believe that the problematic theoretical 
density dependence of conductivity for low N Q is associ- 
ated with the fact that, in the dilute limit, the response of 
the system caused by the electric field can not be treated 
as linear in |E|. It is well known that linear response 
theory is valid only for e|E|a < jk F with a =< r jt 
as the average displacement induced by the electric field 
(r is the relaxation time). In the dilute limit, k F is very 
small and the requirement e|E|a < jkp is difficult to sat- 
isfy experimentally. Hence, our theory, as well as other 
linear response theories, can not be employed to describe 
the transport behavior of carriers in real systems with 



IV. CONCLUSIONS 



We have formulated a kinetic equation to investigate 
dc transport in graphene considering long-range electron- 
impurity scatterings. In this study, we included the role 
of interband correlations. We found that the conductiv- 
ity contains two terms: one term is inversely proportional 
to the impurity density while the other one is linear in 
Ni. This results in a minimum in the density depen- 
dence of conductivity for N e equal to a critical density, 
N c . For N c > N c the conductivity varies almost linearly 
with the electron density while, for N c < N c , a is approx- 
imately inversely proportional to N c . We also discussed 
the effects of various scattering potentials on the conduc- 
tivity minimum. Using typical experimental parameters, 
we found that for the RPA-screened Coulomb potential, 
cr min w 5.1 e 2 /h and N c w 0.327V,. 



Furthermore, the asymptotic dependence of conductivity 
on electron density as N c moves away from the critical 
value is given in the TF and RPA cases as: 



iTF 
\T=0 



1 RPA 




20.07VC/7V 
0.257V/A C 



■ONc/Nt 
Q.8Ni/N e 



TV 



>NJ F 

< nJ f 



N c > iV c RPA 
N r < 7V RPA 



(42) 



(43) 
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N c /Ni 


o- m in/(e 2 /ft) 


SR 7iV<5/32 « 0.24 


2^6 « 4.9 


TF v / F(2r s )G(2r 3 )/2 


2[VF(2r s )/G(2r a )] 


RPA /F[4r s /(2 - 7ir s )]C?[4r s /(2 - 


- vrr s )]/2 2{ 1 /J ? [4r s /(2 - 7rr a )]/G[4r 3 /(2 - 7rr s )]} 



TABLE I: Critical electron densities and minimum conductivities for the various screened Coulomb potentials considered 
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